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Abstract 

A short review of experimental and theoretical results on the large angle cross sections 
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"77 — y two mesons" and the form factors 7*7 — > P = {11, 77, 77'} is given. 



1. Introduction 

The general approach to calculations of hard exclusive processes in QCD was developed in 
[HE]. In particular, the general formula for the leading power term of any hadron form factor 



7* — > HiH 2 has the form pQ 

/ . \ |Ai+A2| + (2n min — 3) 

(Pi, Ai; P2, s 2 , A 2 |J A |0) = C 12 (l/\/q 2 ) , (1) 

> , 

^ where : n m \ n is the minimal number of elementary constituents in a given hadron, n m i n = 2 for 
v^q mesons and n m i n = 3 for baryons ; 51,2 and Ai^ are the hadron spins and helicities, the current 
O I helicity A = Ai — A2 = 0, ±1 ; the coefficient Cyi is expressed through the integral over the wave 
CN " functions of both hadrons. 

Q\ \ It is seen that the behavior is independent of hadron spins, but depends essentially on their 
helicities, and the QCD helicity selection rules are clearly seen: the largest form factors occur 
only for Ai = A2 = mesons and Ai = — A2 = ±1/2 baryons of any spins. 

^ 1 The QCD logarithmic loop corrections to (1) were first calculated in [2] (see also [3]jl][5] for 
■ more details and [6] for a review). 

2. 77 — >• MM large angle scattering 

The QCD predictions for the leading terms of the large angle scattering cross sections 77 — >• 
two mesons were considered in [7] [8] (see also [9] for the one-loop corrections) 





Fig.l Two typical lowest order Fcynman diagrams for the leading term hard QCD contributions 
to 77 — > MM , the broken line is the hard gluon exchange. 
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The expressions for the cross sections look as (the example in (2) is given for 77 — > K + K ) : 



da(77 -»■ M f M) 
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where: x s + x u = 1 , e u = 2/3, e s = = —1/3 , fp are the couplings : f n ~ 132 MeV , ~ 
162 MeV, 4>p(x) is the leading twist pseudoscalar meson wave function (= distribution ampli- 
tude), "x" is the meson momentum fraction carried by quark inside the meson. 




Fig. 2 Three different models for the leading twist pion wave function cji-^ix). 
Red line - asymptotic wave function cj)% By (x) = 6x(l — x) . 
Blue line - CZ wave function (at the low scale normalization point £t D ~ lGeV) : <j>° z (x, ^t ) = 30x(l — x)(2x — l) 2 |10j . 

Black line - flat wave function <f> n (x,fi ~ 1 GeV) = 1. 

Cross sections for charged mesons: 77 —> 7r + 7r~, K + K~ behave as: 

Ar(77 -» 7T+7T-) f* 

dcos8 ~ jy 6 sin 4 # ' 1 j 

and the angular distribution ~ 1/ sin 4 9 is only weakly dependent of the meson wave function 
form. But the absolute values of cross sections depend strongly on the form of 4>m(x) and are 
much larger for the wide wave functions. 

For neutral mesons: 77 — > 7r°7r°, K$Ks, 7T°r], r\r\ the coefficient of the formally leading term 
~ 1/W 6 is very small, so that at present energies W < A GeV such amplitudes are dominated by 
the first power correction in the amplitude and the energy behavior is much steeper: 
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(4) 



while, unlike (3), the angular dependence and the overall coefficient in (4) are not predicted 
(at present) in a model independent way. 

As the alternative approach to description of 77 — > MM processes, the "handbag model" was 
used in [11]. The main dynamical assumption of the "handbag model" [11] is that at present 
energies W < 4 GeV all 77 — > MM amplitudes are still dominated by soft non-leading terms. 




Fig. 3a The overall picture of the "standard handbag" contribution 1111 . 



Fig. 3b The standard lowest order Feynman diagram for the QCD light cone sum rule 1121 . 



For all mesons, both charged and neutral, "the standard handbag" contribution (Fig. 3) gives: 
rfo"(77 — > MM)/d cos 9 ~ const/H^ 10 [12]. This angular behavior ~ const disagrees with all data 
~ 1/ sin 4 9, and the energy behavior disagrees with the data ~ 1 /W 6 for charged mesons. 
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Fig. 4 The typical additional Feynman diagram 
for the "extended handbag model" which includes contributions 
from 3-particle wave functions (the curly line is the near mass 
shell non-pcrturbative gluon). 



I expect that, in distinction with the standard contribution of Fig.3, the additional soft con- 
tributions like those in Fig.4 will give : 
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in better qualitative agreement with data for neutral mesons. Unfortunately, such contributions 
are not yet calculated at present (and it well may be that they are too small in absolute values). 

Now, about a comparison with the data. The Belle results for 77 — > tt + tt~ and 77 — > K + K~ 
are presented in Fig. 5. 
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Fig. 5 a,b) Cross sections <j integrated over the angular region | cos^l < 0.6 , together with the ~ (l/W)® dependence line [14] . 
c) the cross section ratio iJ cxp = o- (K J '~ K~)/o- (tt + n~ ) ~ 0.9 1141 . Compare -R exp — 0.9 with the naive prediction R = {fit/fir) 4 — 
2.3. 1 



1 It was "obtained" in [TT| that the angular behavior of the standard handbag contribution from Fig. 3a is 
da/dcos9 ~ 1/sin 9. Really, this "result" is completely model dependent. The reason is that a number of special 
approximate relations were used in at intermediate steps. All these relations are valid, at best, for the leading 
term only. But it turned out finally that their would be leading term gives zero contribution to the amplitude, 
and the whole answer is due to next corrections, which were not under control in Their result ~ 1/ sin 9 
is completely due to especially (and arbitrary) chosen form of the next to leading correction, while ignoring all 
others next to leading corrections of the same order of smallncss. Therefore, there is no really model independent 
prediction of the angular dependence in So, it is not surprising that the explicit calculation in [12| gives 
different angular dependence. 

2 The reason for i? oxp <C (fx/ fw) 4 is that the leading twist pseudoscalar meson wave function <fip(x) becomes 
narrower when the lighter u or d quarks are replaced with the heavier s quarks, and this opposite effect compensates 
those from fa/ fir > 1, see [8]. 



3 



It is seen that they are compatible with the leading term QCD predictions a(ir + 7T ) ~ 
a(K + K~) ~ 1/W 6 , and disagree with the handbag model predictions of much steeper behavior 
~ 1/W W . 
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Fig. 6 (a) The total cross section cr (77 ■ 
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KgKg) integrated over the angular region | cos6| < 0.6 1151 . 



Here n is the TV-dependence a (W) ~ 1/W™; (b) The ratio a (K s K s )/a Q (K +K~) versus W. 

The dotted line DKV = Diehl-Kroll-Vogt is the handbag model prediction in the SU(3) -flavor symmetry limit 1111 ; the dashed 
BL line is the Brodsky-Lcpage [7] prediction for the kaon wave function close to a sy(a;), the dashed-dotted BC line is the Benayoun- 
Chernyak [8] prediction for the kaon wave function like <f> cz ,{x) (both are the leading term QCD predictions (2) for large energies W). 



The angular behavior measured by Belle [H] for and K + K~ is ~ 1/ sin 4 #, also in 

agreement with QCD and in disagreement with the standard handbag model [12]. As for the 
absolute values of cross sections, the values predicted from (2) are much smaller than data for 
the pion (kaon) wave functions close to (pasyix), while predictions from (2) for the wide n and K 
wave functions like 4> C7 ,{x) are in a reasonable agreement with data (see [12] for more details). 

Now, let us compare with the Belle results for the neutral mesons. The results for the cross 
section 77 — > K$Ks are published in [15], see Fig. 6. 

It is seen from Fig.6 that in the energy range 2.5 < W < AGeV the energy behavior ~ 1/M^ 10 
in this neutral channel is much steeper in comparison with ~ 1/W 6 in the charged channel. This 
agrees with qualitative expectations from QCD that because the coefficient of the formally leading 
at sufficiently large W is very small, the first non-leading term dominates the i^i^s-amplitude 
at present energies W < 4 GeV . Let us recall that the handbag model predicts the dominance 
of non-leading terms (and so the energy behavior ~ 1/W 10 ) for all mesons, both charged and 
neutral. 

As for the angular distribution, the data are sufficiently well described by ~ 1/ sin 9 |15j . 
Let us recall once more that "the standard handbag model" (Fig. 3) predicts the flat angular 
distribution ~ const [12] also for all mesons, but the qualitative expectation is that "the extended 
handbag model" (Fig. 4) will give ~ 1/ sin 4 #, see (5). 

Finally, about the ratio R = a (KsKs) /o~q(K + K~) , see Fig.6b. In the SU(3) flavor symmetry 
limit the standard (and extended) handbag model predicts R = 0.08 [IT] . It is seen from Fig. 6b 
that this ratio decreases rapidly with energy and becomes smaller than ~ 0.08 at W > 2.7 GeV, 
in disagreement with the handbag model. This is because the energy dependence of o~q{K + K _ ) ~ 
1/W 6 disagrees with the handbag model prediction ~ l/W 10 . 

The QCD prediction is that at sufficiently large W, when the parametrically leading but 
having a small coefficient term will become dominant in the K$Ks amplitude, this ratio will 
become constant (see BL and BC lines in Fig. 6b). It is seen from Fig. 6b that the ratio R is 
already close to the leading term QCD predictions for KsK$ at W ~ 4 GeV . 

The qualitative situation with other neutral modes, 77 — > it ir°, ir°ri, rjrj, rj't]', etc., is similar 
to those of 77 — > KgKs- Recently, there appeared new data from the Belle Collaboration on cross 
sections 77 — > 7i ir and 77 — > 7r rj [TB] [T7], see Figs. 7,8 and Table 1. The QCD predictions for this 
range of energies are: ct(7t + 7t~) ~ 1/W 6 , a(7r°iv ) ~ a{TX r]) ~ 1/W 10 , R = o{tx°ix ) / o(ix + ix~) ~ 
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1/W A . The handbag model prediction is: R = a(7r 7r )/cr(7r + 7r~) = 1/2 . As for the cross section 
a(77 — > TT r]), it behaves "normally", ~ 1/W W , similarly to cr(77 —> KsKs), see Fig. 8a. But 
as for a (77 — > 7r°7r°), it behaves "abnormally" , see Fig. 7. This last behavior agrees neither with 
QCD, nor with the handbag model. 

3. Conclusions on the large angle cross sections 77 — >■ MM 

1) The leading term QCD predictions da/d cos 9 ~ l/fW 6 sin 4 6) for charged mesons 7r + 7r~, K + K~ 
agree sufficiently well with data both in energy and angular dependence at energies W > 3 GeV . 
The absolute values of cross sections are in a reasonable agreement with data only for the wide 

it (K) wave functions, like <\f* K {x) (see [8j [12] for more details). The asymptotic wave func- 
tions (fr-x^ix) — 4> asy {x) predict much smaller cross sections. The handbag model predictions for 
charged mesons disagree with data in energy dependence. 

2) For neutral mesons, the QCD leading terms have much smaller overall coefficients, so that 
the non-leading terms are expected to dominate at present energies W < 4 GeV, and the energy 
dependence is steeper: a(M°M°) ~ 1/W 10 . This agrees with data on a(KsK s ) and a(ir°r]), 
while cr(7r°7r°) behaves "abnormally" (may be due to contamination of data with the pure QED 
- background). 

3) Predictions of the "standard handbag model" disagree with data either in energy and/or 
angular dependence, or in absolute values. However, it is not excluded that adding soft contribu- 
tions from the 3-particle wave functions in the "extended handbag model" (see Fig.4) can help 
to describe cross sections of neutral mesons at intermediate energies 2.5 GeV < W < 4. Unfor- 
tunately, such contributions are not yet calculated at present (and it well may be that they are 
too small; besides, one has to remember that there are also power corrections due to the higher 
twist wave function components in the diagrams in Fig.l ). 

Table 1 : The value of "n" in otot ~ (1/W) n in various reactions fitted in the W and | cos 8\ ranges indicated. 



Process 


n - experiment 


W range (GeV) 


cos 9 


Ref. 


n - QCD 


n - handbag 


7T + 7T~ 


7.9 ± 0.4 stat ± 1.5 syst 


3.0-4.1 


< 0.6 


m 


~ 6 


~ 10 


K + K~ 


7.3 ± 0.3 s tat i 1.5 sys t 


3.0-4.1 


< 0.6 


114 


~ 6 


~ 10 




10.5 ± 0.6 stat ± 0.5 syst 


2.4 - 4.0 


< 0.6 


115 


~ 10 


~ 10 


r]7i° 


10.5±1.2 stat ±0.5 syst 


3.1 - 4.1 


< 0.8 


117 


~ 10 


~ 10 


7T 7r° 


8.0 ± 0.5 stat ± 0.4 syst ? 


3.1 - 4.1 


< 0.8 


116 


~ 10 


~ 10 
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Fig. 7 (a) Cross sections cr (77 — > tt°tt°) 
and 0-0(77 -> 7T+7T-) for | cos6»*| < 0.6 [161 [14| . 
(b) Their ratio. 
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8 a) W - dependence of cross sections 77 — > tt°tt° and 77 — > tt°T] (| cos6*| < 0.< 

The curve is the fit: a(Tr°~q) ~ W~ n , n = (10.5 ± 1.2 ± 0.5) [TBIfTT] . 
b) W - dependence of the cross section ratio a(rjir°) /a(n°iT ) (| cos9*| < 0.8). 



4. 7*7 — ► P = "t 71 " j *7 ? »7'} form factors i^yp(Q 2 ) 



As was first obtained in j2] on the example of the pion form factor F n (Q 2 ) (see [3] for more 
details), the contributions from short and large distances factorize in F n (Q 2 ) at large Q 2 , and the 
logarithmic evolution of Q 2 F n (Q 2 ) is determined by renormalization factors of operators with the 
same anomalous dimensions 7„ as in the deep inelastic scattering. So, the strict non-perturbative 
QCD prediction for F^^Q 2 ) in the formal limit Q 2 — > oo looks as [2] (b = 11— 2n//3, /i ~ 1 GeF) : 



This corresponds to the pion wave function (j) n (x,fi — > oo) evolving to its universal asymptotic 

f0mi 7r (x,/i>lGeV)->6x(l-x)^l + o(^^)^j , 

independently of its form (f) w (x, /x G ~ 1 GeV) at low energy. As it is seen, the logarithmic evolution 
with increasing scale is very slow. 

As for the form factor F J7T (Q 2 ) = F~ /7T (Q 2 = —qf, q\ = 0), the QCD prediction [f| for its 
asymptotic behavior in the formal limit Q 2 — > oo looks as (see e.g. [5]) : 

J dze^ z {-K(p)\T{J^z)J u (Q)}\®) = (i^uxaq^) F J7T (Q 2 ), (6) 

For the r] and rf mesons, the form factors F ir] and F 1V ' look similarly to F ln . For instance, a 
simplified description of \rj), \rf) states in the quark flavor basis looks as follows [i~9] : 

\tt°) ->■ |(mt - dd)/V2), \n) ->• |(«« + dd)/V2), \s)-+\ss), 

\rj) = cos0 |n) — sin0 |s) , |r/) = sin |n) + cos |s) . (7) 
U~132MeV, U^U, fs^l-3U, 0-38°. 

F„(C?)= ^-® f * [\ ^ X >»~ Q h o] F in {Q 2 )= ^X e2d)U fd* UX ^~ Q) Io, 
Q Jo x Q Jo x 

F, S (Q 2 ) = ^ j 1 dx 0a(x '^ g) I ; Jo = ( 1 + 0(«.) + 0(l/g 2 ) j , g 2 » 1 Ge\/ 2 . (8) 

F 7?? (g 2 )= |ros</>F 7n (g 2 ) -sin0F 7S (g 2 )j, F 7r/ (g 2 )= ^sin0F 7n (g 2 ) + cos0F 7S (g 2 

Predictions for _F 77r (g 2 ) were given in a large number of theoretical papers, using many different 
models for the leading twist pion wave function ^(x, /S). The previous data for F 17T (Q 2 )^20\ [2Tj 
covered the space-like region < Q 2 < 8GeV 2 . The recent data from BaBar [231 121] extended 
this one to Q 2 < 40 GeV 2 , see Fig. 9. It is seen that Q 2 F 1 p{Q 2 ) exceeds its asymptotic value a/2 
(the dashed line in Fig.9) at Q 2 > lOGe^ 2 . Because the loop and leading power corrections are 
negative here, this shows that the leading twist pion wave function <p n (x,fi) is considerably wider 



3 Really, unlike F n (Q 2 ), the leading asymptotic behavior of F 1V (Q 2 ) can be directly obtained from the standard 
Wilson operator expansion of J^{z)J y {Q) ™» J2 n C n (z)O n (Q) in (6) [IS], as in calculations of the deep inelastic 
scattering. The only difference is that the forward matrix elements (p\O n \p) are taken in the deep inelastic 
scattering, while these are (p|O„|0) in the case 7*7 — > ir. 
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than 4>asy(x), while most theoretical models predicted cj) w (x,(j,) ~ <pasy(x). The red curve in Fig.9 
shows that, with power corrections of reasonable size, the wide leading twist pion wave function 
4> cz (x, fi) obtained in [10] using the standard QCD sum rules, see Fig. 2, is not in contradiction 
with data. 

After the new BaBar data on F 77T (Q 2 ) [23] appeared, it was proposed in [28] that the large 
value of Q 2 F- in (Q 2 ) is due to the flat pion wave function, (j) n (x,fi ~ 1 GeV) ~ 1, as Q 2 F ln (Q 2 ) 
grows ~ ln(Q 2 /M 2 ) in this case and, taking M ~ m p by hand, such a behavior fits well then these 
BaBar data. But in this case, because the wave functions of \ir), \n) and \s) are qualitatively 
similar, see (7-8), the form factors q 2 F~ fri (q 2 ) and q 2 F~ /r) i(q 2 ) will also grow the same way, ~ 
\n(q 2 /m 2 ) at q 2 3> lGeV 2 . These two form factors have been measured recently by the BaBar 
Collaboration [22] at q 2 = 112 GeV 2 . It is seen from Table 2 and Fig. 10 that with 4> n (x) ~ <f> a { x ) — 
1 these two form factors will be too large. Besides, such flat wave functions will contradict the 
data on <r(e + e~ — > VP), see Fig. 11 . 

Tabic 2. The values of form factors \q 2 F 7 p(q 2 )\ (in GeV) at q 2 = 112 GeV 2 for various meson wave functions 



Wave functions 


\q 2 F r ^(q 2 )\ 


\Q 2 F rn (q 2 )\ 


\q 2 F rv '(q 2 )\ 


Ref. 


<j>n(x) ~ (f> s (x) ~ <t>asy{x) = 6x(l - x) 


0.14 


0.13 


0.21 




4> n (x) ~ 4> s (x) ~ (f) cz {x) 


0.22 


0.21 


0.33 




0„(x) ~ cz (a;); <fi s (x) ~ asy (a;) 


0.22 


0.24 


0.29 




<j>n(x) ~ (f> s (x) ~ 1 


0.32 


0.31 


0.49 




experiment 




0.23 ±0.03 


0.25 ±0.02 


[21 



4 A number of papers with predictions for F 17T {Q 2 ) has been published, based on the model pion wave function 
(f>® MS (x) (BMS=Bakulev-Mikhailov-Stefanis) obtained from the "improved QCD sum rules" with non-local con- 
densates (see the last paper [25] and references therein) . This approach has been criticized in [13] , as it is based on 
arbitrary strong dynamical assumptions which, as was shown in |13j . don't pass the direct QCD check. Moreover, 
within this approach, one has to introduce a number of arbitrary model functions for various non-local vacuum 
condensates (see e.g. [26]) and, in general, the results for (j) n (x) depend heavily on the model forms chosen for these 
functions (compare e.g. the results for ff (x) from [37] and US]). Finally, the model pion wave function 4>® MS (x, /i) 
obtained within this approach predicted the value of F^(Q 2 ) only slightly above those for the asymptotic wave 
function ^(x) = </> asy (:r), see [25], and well below the recent BaBar data [23] [24]. Besides, it is claimed in [25] 
that the data [20] [21] [23] [24] are incompatible with (f) n (x,fj) — cf) cz (x,fj,) (and are even in contradiction with the 
QCD factorization for any pion wave function with the end point behavior ~ x(l — x) at x — > 0, 1 ). As it is seen 
from Fig.9 (red curve), this is not so. 
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Fig.9 The form factor <3> : 



□ CELLO 
A CLEO 
• BABAR 




Theory : 

a) the logarithmic loop corrections are calculated (in part) at the NNLO [29 30 ; 

b) only the part of the total power correction ~ 1/Q 2 is calculated at present in 1311 : 

5<t> 4 ~ -V2U(0.6GeV 2 )/Q 2 , 

originating from the 2 - and 3 - particle asymptotic pion wave functions of twist 4 1321 . 
It well may be that it is not even the main part of the total ~ 1/Q 2 correction, because: 
i) the deviation of these 2 - and 3 - particle twist 4 pion wave functions from their asymptotic 
forms (i.e. the admixture of higher "partial waves") was estimated in |32] only for the first 
non-leading "partial wave", while second (and higher) "partial waves" neglected in 1321 can be 
important, ii) there are also contributions ~ 1/Q 2 from the 4 -particle wave functions of 
twist 4, iii) moreover, the twist expansion breaks down at this level, so that the higher 
twist > 6 terms also give contributions ~ 1/Q 2 ; 
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c) the power correction ~ 1/Q 4 is unknown. 

Black line: 0^(z) = </> asy (x), * ~ v^jV [ 0.77 - (0.6 GeV 2 IQ 2 ) ] [29 ] [30 ] [31] 
Blue line : ^(x) = cf> C z(x), <S> ~ %/2/V [ 1.18 - (0.6 GeV 2 /Q 2 ) ] [30 ] [31 ] 
Red line (the example with additional mild power corrections) : 
<MaO = <£cz(£c), * Oi V*f* [ 1.18- (1.5 GeV 2 /Q 2 ) - (1.2 GeV 2 /Q 2 )* 

Experiment : 




10 

q 2 (GeV 2 ) 



q (GeV) 



Fig. 10 Full points: \q 2 F ir] (q 2 )\ (left) and \q 2 F jr) i (q 2 )\ (right) transition form factors at q 2 = \\2GeV 2 
\q 2 F lv {q 2 )\ = (0.229 ± 0.030 ± 0.008) GeV, \q 2 F irj ,{q 2 )\ = (0.251 ± 0.019 ± 0.008) GeV. 

White points: previous CLEO data at 2 GeV 2 < (Q 2 = -q 2 ) < 20 GeV 2 [3T]. 

Black lines: the form factors \q 2 F JT) (q 2 )\, \q 2 F J7] > (q 2 )\ for the fiat pseudoscalar wave function (j>p( x ) — 1 



5. Conclusions on the form factors F lP (Q 2 ), P = rj, rj'} 
and the leading twist wave functions <fip(x) of pseudoscalar mesons 

The flat leading twist pseudoscalar wave function <pp(x) ~ 1 : 

a) predicts the form factors F lv (q 2 ) and F lv /(q 2 ) at q 2 = 112 GeV 2 considerably larger than the 
BaBar results; 

b) predicts the parametrical behavior of cross sections a(e + e~ — > VP) at large s as: cr(e + e~ — > 
VP) ~ 1/s 2 , in contradiction with the data a(e + e~ — > VP) ~ 1/s 4 in the interval ~ 8 GeV 2 < 
s < lllGeV 2 . 

The asymptotic leading twist pseudoscalar wave function 4>p(x) ~ 6x(l — x) : 

a) predicts the form factors F 77r o(Q 2 ), F 1V (Q 2 ) considerably smaller than data; 

b) predicts branchings of charmonium decays: 3 Pq, 3 P2 — > tt + it^, K + K~, the pion and kaon 
electromagnetic form factors F^^io 2 ) at q 2 = 10 — 15 GeV 2 much smaller than data, etc. 
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The CZ leading twist pion wave function 4>^ z {x, \x ~ 1 GeV) = 30x(l — x)(2x — l) 2 
leads to predictions which, it seems, are not in contradiction with all data available. 



Fig. 11 Solid lines correspond to 1/s 4 dependence 
and dashed ones represent 1/s 3 . 

a) a(e + e~ — > 0ry) b) cr(e + e~ — > 4>t]') 
c) cr(e + e~ — > prj) d) cr(e + e" — > prj') 

The measured cross sections: 
at v^s ~ 2.5, 2.75 GeV by BaBar [33], 
at y/s = 3.67 GeV by CLEO J3J, 
at yfs = 10.58 GeV by BaBar [55] and Belle [35] 
for various processes. 

BaBar measurements are represented by squares . 

QCD predictions : 

cr(e+e- -> VP) ~ 1/s 4 [T| ( see (1), |A V | = 1 in this case), 
up to a possible additional logarithmically growing factor [6], 
for the pseudoscalar wave function <f>p(x) with the suppressed 
end point behavior, like ~ x(l — x) at x — > 0, 1. 

<x(e+e — — > VP) ~ 1/s 2 for the flat pseudoscalar 
wave function <j>p(x) ~ 1. 

The form factors 7* — > 1/P are highly sensitive to the end point behavior of the leading twist pseudoscalar wave function <f>p(x), 
as they contain the factor / ~ Jg dx </>p(x)/x 2 , <5 = 0(^i„/Q 2 ) [6]. So, I ~ ln(Q 2 /^t 2 ) for <j>p(x) ~ x(l — x), while it will be 
parametrically larger at <j>p( x ) ~ 1 : / ~ (Q 2 //i 2 ). The data are in a reasonable agreement (with a logarithmic accuracy) with the 
a ~ 1/s 4 dependence corresponding the end point behavior <pp(x) ~ x(l — 2), and are in contradiction with the behavior a ~ 1/s 2 , 
corresponding to <pp(x) ~ 1 at i -> 0, 1. 

Note added 

After this talk has been given, there appeared the paper [37j with updated "predictions" of 
the handbag model for the 77 — > MM cross sections. In comparison with the previous paper [11] , 
the main new element in [37] is that the sizeable soft non- valence form factor R^„Js) is used 

1 1 r-i MM K > 

now, in addition to the soft valence one, BZj (s). □ Both functions, (s) and (s), are 

IVl Ivl 1Y1 ±\{l 1\± 

parameterized then in arbitrary forms, with a large number of free parameters which are fitted 
to the data. It seems, the authors consider that nearly nothing is known about QCD, so that 
it is possible to proceed in such an arbitrary way. 

As for the standard soft valence contributions to the cross sections and the soft valence form 
factors K^j M (s), these were estimated in [32] with a help of the QCD light cone sum rules (see 
section 2 above) and were found much smaller (and with the behaviour ~ 1/s 2 ) than the values 
fitted in [TT] (and in [37]). In addition, we would like to comment here in short on the non- valence 
contributions. 

5 As for the "prediction" of the angular behavior sda/d cos 9 ~ \Rjj m \ 2 /sin 4 9 in [11], see footnote 1 and [12] . 

6 The form factors R^{s) and R s 2jr (s) used in [37] are: i? 2 " 7r (s) = R v 2?r (s) + i?^(s), R s 2 „(s) = R^l(s). 
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Two types of non- valence contributions are presented in Fig. 12. 



Fig. 12a The leading power non- valence one- loop correction. 
Fig. 12b The contribution to the soft handbag form factor _R2X^(s). 
The solid and dashed lines represent quarks and gluons. 

It is worth noting that both non- valence contributions in Fig. 12 are SU (3)-flavor singlets in the 
5'(7(3)-symmetry limit. So, they contribute equally to the amplitudes 7r + 7r~, 7r°7r°, K + K~ , K°K° 
and r7 8 ?78, and don't contribute to r] 8 7r°. 

The diagrams in Fig. 12a constitute a small subset of all one- loop corrections to the leading 
power contributions from the Born diagrams like those shown in Fig.l . We only note here that 
if these leading power one-loop non-valence corrections to the Born contributions were really 
significant, this will contradict then the data [15] on K$Ks, see Fig.6. 

As for the soft non- valence handbag form factor R^ M (s), it seems sufficient to say that it orig- 
inates first from the Fig. 12b two-loop (non-logarithmic) correction, so that : R^ M (s) / R^ M (s) = 

o(^(a s (s)/7r) 2 ^j = O(10~ 2 ) . Clearly, so small non-valence contribution will not help. 

I am grateful to V.P. Druzhinin for explaining me details of various experimental results and 
G. Duplancic for providing me with the additional details of calculations performed in [9]. 

This work is supported in part by the RFBR grant 07-02-00361-a. 
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